Quantum state filtering and discrimination between sets of Boolean functions 
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In quantum state filtering one wants to determine whether an unknown quantum state, which 
is chosen from a known set of states, • • • ji/'iv)}, is either a specific state, say or one of 

the remaining states, {IV'2), • ■ • \4'n)}- We present the optimal solution to this problem, in terms 
of generalized measurements, for the case that the filtering is required to be unambiguous. As an 
application, we propose an efficient, probabilistic quantum algorithm for distinguishing between sets 
of Boolean functions, which is a generalization of the Deutsch-Jozsa algorithm. 

PACS numbers: 03.67-a, 03.65.Ta, 42.50.-p 



Optimal discrimination among quantum states plays 
a central role in quantum information theory. Interest 
in this problem was prompted by the suggestion to use 
nonorthogonal quantum states for communication in cer- 
tain secure quantum cryptographic protocols, most no- 
tably in the one based on the two-state procedure as de- 
veloped by Bennett [J. The reason why until recently 
the area has shown relatively slow progress within the 
rapidly evolving field of quantum information is that it 
poses quite formidable mathematical challenges. Except 
for a handful of very special cases, no general exact so- 
lution has been available involving more than two arbi- 
trary states. In this paper we present an exact solution 
to an optimum measurement problem involving an arbi- 
trary number of quantum states, with no restriction on 
the states. The resulting method has the potential for 
widespread applications in quantum information process- 
ing. In particular, it lends itself quite naturally to a quan- 
tum generalization of probabilistic classical algorithms. 
Whenever it is possible to find a one-to-one mapping of 
classical alternatives onto quantum states, our method 
can discriminate among these quantum alternatives in a 
single step with optimum success probability. 

We illustrate the strength of the method on the exam- 
ple of a probabilistic quantum algorithm to discriminate 
between sets of Boolean functions. A Boolean function on 
n bits is one that returns either or 1 as output for every 
possible value of the input x, where < a: < 2" — 1. The 
function is uniform (or constant) if it returns the same 
output on all of its arguments, i.e. either all O's or all 
I's; it is balanced (or even) if it returns O's on half of its 
arguments and I's on the other half; and it is biased if it 
returns O's on mp of its arguments and I's on the remain- 
ing TOi = 2" — mo arguments (mp ^ mi 7^ or 2" — 1). 
Classically, if one is given an unknown function and told 
that it is either balanced or uniform, one needs 2("-i)-hl 
measurements to decide which. Deutsch and Jozsa 
developed a quantum algorithm that can accomplish this 
task in one step. To discriminate a biased Boolean func- 
tion from an unknown balanced one, 2^""-'^) -I- mi + 1 



measurements are needed classically, where, without loss 
of generality, we have assumed that mi < mg. Here we 
propose a probabilistic quantum algorithm that can un- 
ambiguously discriminate a known biased Boolean func- 
tion from a given set of balanced ones in one step. 

The method is based on the optimum unambiguous 
quantum state filtering scenario which, in turn, is a spe- 
cial case of the following more general problem. We know 
that a given system is prepared in one of N known non- 
orthogonal quantum states, but we do not know which 
one. We want to assign the state of this system to one 
or the other of two complementary subsets of the set of 
the N given states where one subset has M elements and 
the other has N - M (M < N/2). Since the subsets 
are not mutually orthogonal, the assignment can not be 
done with a 100% probability of success. For the case 
that the assignment is required to be unambiguous, at 
the expense of allowing inconclusive results to occur the 
probability of which is minimized, the problem has re- 
cently been solved for = 3 3]. For the case that the 
assignment is to be performed with minimum error, the 
solution has been found for arbitrary M and N under the 
restriction that the Hilbert space spanned by the states is 
two-dimensional 01 • We refer to either case as quantum 
state filtering when M = 1 and A^ > 3. 

Unambiguous filtering is related to unambiguous quan- 
tum state discrimination: one is given a quantum system, 
prepared in a state that belongs to a known set of non- 
orthogonal states, and one wants to determinCjWithout 
possibility of error, which state the system is in p. Since 
the states are not mutually orthogonal, at first glance the 
problem appears impossible. However, it becomes pos- 
sible if we allow the procedure to fail a certain fraction 
of the time. That is, when we apply the procedure, we 
will either find out what the quantum state of the sys- 
tem is, or we will fail to do so, but we will never make 
an erroneous identification. The optimal method for dis- 
criminating between two states was found in Refs. ^J-^]. 
No general solution is known for more than two states but 
there are special cases that can be solved, and some fea- 
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tures of the general solution have been extracted |9j-[13|- 
Chefles has shown that the states have to be lin- 
early independent for unambiguous state discrimination 
to succeed but orthogonality is not required. 

We begin by presenting the solution to the unam- 
biguous quantum state filtering problem. Suppose we 
are given a quantum system prepared in the state 
which is guaranteed to be a member of the set of N 
non-orthogonal states , . . . ,\iIjn)}, but wc do not 
know which one. We denote by r]i the a priori prob- 
ability that the system was prepared in the state 
We want to find a procedure that will unambiguously 
assign the state of the quantum system to one or the 
other of two complementary subsets of the set of the N 
given non-orthogonal quantum states, either {I'f/'i)} or 
{1^2), ■■ ■ \^n)}- Quantum measurement theory tells us 
that non-orthogonal states cannot be discriminated per- 
fectly. If we are given we will have some probability 
Pi to correctly assign it to one of the subsets and, cor- 
respondingly, some failure probability, qi = 1 — pi, to 
obtain an inconclusive answer. The average probabilities 
of success and of failure are P = X^t^i ViPij ^^id 

N 

Q = ^Viqi, (1) 

respectively. Our objective is to find the set of {qi} that 
minimizes the probability of failure, Q. 

It is easy to see that a standard quantum measure- 
ment (SQM, a von Neumann projective measurement) 
can achieve error-free filtering. If we project on cither of 
the two sets (state selective measurement, first strategy) 
a "no click" will indicate that we were given a state from 
the other set, assuming perfect detectors. A somewhat 
better approach is to project on a direction that is per- 
pendicular to one of the sets (nonselective measurement, 
second strategy). Now, a detector "click" indicates that 
we were given a state from the other set and perfect de- 
tectors are not required. For example, if we measure the 
operator F^^^ = / — \ipi){'tpi\, then a click (correspond- 
ing to the eigenvalue 1) shows that the vector is not 
and the measurement has succeeded. If we do not obtain 
a click (eigenvalue 0), then the measurement has failed, 
and we do not know which vector we were given. The 
probability of failure, Q'^q^, is given by 

QSQM =Vi+S, (2) 

where S — X)iL2 ^iKV'ilV'i)!^ is the average overlap be- 
tween the two subsets. 

A second possibility is to split into two com- 

ponents, Itjji) = + lipi)- Here {-ip^) is orthogo- 

nal to the subspace, 7^2, that is spanned by the vec- 
tors 1-02), ... li/^Af), and IV'i) lies in H2- Their normahzed 
versions are {ipi) = {ipi) /\\^i\\ and I'tpl) = |i/'i)/ll V'l II > 



respectively, where the norm is defined in the usual 
way, IIV'lP = {ip\ip)- We then introduce the operator 
= - - which has 

eigenvalues 1, 0, and —1. If we measure F^^^ and obtain 
1, then the vector was if we obtain —1, then the 

vector was in the set {\ip2), ■ ■ - IiPn)}, and if we obtain 

0, the procedure failed. In this case the probability of 

(2) 

failure, QgQ}^J^ is given by 

QfQM = m\\4r + -A-. (3) 
llV-ilr 

Which of these two particular strategies is better is de- 
termined by which of these two failure probabilities is 

smaller. In particular, Q^^q^ > Q^sqm '7i|IV'ilP > 
and vice versa. 

Now, the question arises: Is this the best we can do? 
The answer is that under certain conditions a generalized 
measurement based on positive-operator valued measures 
(POVM, [13) can do better in an intermediate range of 
parameters, and can achieve a higher probability of suc- 
cess than a standard von Neumann measurement. The 
POVM can be implemented by a unitary evolution on 
a larger space and a selective measurement. The larger 
space consists of two orthogonal subspaces, the original 
system space and a failure space. The unitary evolution 
transforms the input sets into orthogonal sets in the orig- 
inal system space and maps them onto the same vector in 
the failure space. A click in the detector measuring along 
this vector corresponds to failure of the procedure, since 
all inputs are mapped onto the same output. A no-click 
corresponds to success since now the non-orthogonal in- 
put sets are transformed into orthogonal output sets in 
the system space. The one-dimensionality of the failure 
space follows from the requirement that the filtering is 
optimum. Namely, suppose that j-i/'i) is mapped onto 
some vector in the failure space and the inputs from the 
other set are mapped onto vectors that have components 
perpendicular to this vector. Then a single von Neumann 
measurement along the orthogonal direction could iden- 
tify the input as being from the second set, i.e. further fil- 
tering would be possible, lowering the failure probability 
and the original filtering could not have been optimum. 

In particular, let Tig be the £)-dimensional system 
space spanned by the vectors {|V'i), ■ • ■ I'/'w)} where, ob- 
viously, D < N. We now embed this space in a space 
of -D -|- 1 dimensions, Hs+a = ® Ha, where Ha 
is a one-dimensional auxiliary Hilbert space, the failure 
space or ancilla. The basis in this space is denoted by 

[(/i^"^)). Thus, the unitary evolution on Hs+a is specified 

( s) 

by the requirement that for any input state \ipi){= IV'l )) 
{i — 1, . . . , N) the final state has the structure 

\^.)out = U\^,) = VP^IVP^) + Vg^e^'^'l^^^)) . (4) 

From unitarity the relation, Pi + qi — 1, follows. Fur- 
thermore, Pi is the probability that the transformation 
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— * succeeds and qi is the probability that 1-0^) is 
mapped onto the state In order to identify pi and 

Qi with the state-specific success and failure probability 
for quantum filtering we have to require that 



(5) 



for z = 2, . . . A^. We now introduce the operator F^'^^ — 
Wilii^'A - (/(^+-^) - Wi){i^'i\ - |(/)(^))((/)(^)|), which has 
eigenvalues 1,0, and —1. If we measure F^^'^ and obtain 
1, then the input was \'4'l)^ if we obtain —1, then the 
input was from the other the set, and if we obtain 0, the 
procedure failed. 

In order to optimize the POVM, we have to determine 
those values of qi in Eq. Q that yield the smallest av- 
erage failure probability Q. Taking the scalar product of 
U\'4)i) and U\'ij)i) in Eq. igj, and using Eq. lO, gives 



(6) 



for i — 2, . . . , iV, and Eq. (Q) can be cast in the form 
Qili) = lili + S/qi- Unitarity of the transformation U 
delivers the necessary condition that qi must lie in the 
range HV'ilP ^ 91 ^ 1- Details of the derivation, along 
with a discussion of the sufficient conditions for the exis- 
tence of U , will be presented in a future publication |l5| . 
Provided that a POVM-solution exists, the minimum of 
Q{qi) is reached for qi = S/rji and is given by 



'.POVM 



(7) 



Thus, the failure probability for optimal unambiguous 
quantum state filtering can be summarized as 




if mU\\\^<s<m, 

if 5 > 7/1 , 

if 5<77iii^lr. 



(8) 



The first line represents the POVM result, Eq. 0, and it 
gives a smaller failure probability, in its range of validity, 
than the von Neumann measurements, Eqs. (|2J) and Eq. 
® , cf. Fig. □ Outside of the POVM range of validity we 
recover the von Neumann results. It should be noted that 
for these results to hold, unlike for unambiguous state 
discrimination, linear independence of all states is not 
required. Instead, the less stringent requirement of the 
linear independence of the sets is sufficient, in agreement 
with the findings in |l(il |. 

We can now apply this result to distinguishing between 
sets of Boolean functions. Let /(x), where < x < 2"— 1, 
be a Boolean function, i.e. f{x) is either or 1. One of the 
sets we want to consider is a set of balanced functions. In 
our example, the second set has only two members, and 
we shall call it Wk- A function is in Wfc if f{x) — for 
< X < [(2'=-l)/2'=]2" and/(j:) = 1 for [(2'''-l)/2'=]2" < 



X < 2" - 
and f{x) 



1, or if f{x) 
= for [(2'= - 



= 1 for < X < [(2*^ - l)/2'=]2" 
l)/2'=]2" < X < 2" - 1. We now 
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FIG. 1: Failure probability, Q, vs. the average overlap, S. 
Dashed line: Q^^qj^^, dotted line: Q^gg^^ solid line: Qpovm- 
For the figure we used the following representative values: 
rj\ — 0.4 and ~ 0.25. For these the optimal Q is given 

(2) w n ^ Q ^ n mc; by Qpgy^j for 0.025 < S < 0.4 



by Q^sQM for < S < 0.025 



and by Q 



(1) 

SQM 



for 0.4 < S. 



wish to distinguish between the given balanced functions 
and functions in W/c, that is, we are given an unknown 
function that is in one of the two sets, and we want to find 
out which set it is in. We note that the two functions in 
Wfc are biased functions, so that this is a special case of 
a more general problem of distinguishing a set of biased 
functions from balanced functions. 

This is by no means the only example the method can 
handle, but it is a particularly simple one and represents 
a generalization of the Deutsch-Jozsa problem Q|. In 
that case one is given an unknown function that is either 
balanced or constant, and one wants to determine which. 
Classically, in the worst case one would have to evaluate 
the function D/2 + 1 times, where we have set D = 2", 
but in the quantum case only one evaluation is necessary. 
The solution makes use of the unitary mapping 



\x)\y)^\x)\y + f{x)), 



(9) 



where the first state, |x), is an n-qubit state, the second 
state, is a single qubit state, and the addition is 
modulo 2. The state |x), where x is an n-digit binary 
number, is a member of the computational basis for n 
qubits, and the state where y is either or 1, is a 
member of the computational basis for a single qubit. 
In solving the Deutsch-Jozsa problem, this mapping is 
employed in the following way 



D-l 



D-l 



J2 kXlO) - |1» - 5](-l)^(^)|x)(|0) - (10) 



x=0 
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This has the effect of mapping Boolean functions to vec- 
tors in the U-dimensional Hilbert space, Hd, and we 
shall do the same. The final qubit is not entangled with 
the remaining n qubits and can be discarded. The vectors 
J2x=o (~^y^^''\^) that are produced by balanced func- 
tions are orthogonal to those produced by constant func- 
tions. This is why the Deutsch-Jozsa problem is easy 
to solve quantum mechanically. In our case, the vectors 
produced by functions in Wk are not orthogonal to those 
produced by balanced functions. However, unambiguous 
quantum state filtering provides an optimum probabilis- 
tic quantum algorithm for the solution of this problem. 

In order to apply the filtering solution, we note that 
both functions in Wk are mapped, up to an overall 
sign, to the same vector in Tijj, which we shall call 
\wk)- The vectors that correspond to balanced func- 
tions are contained in the subspace, Tib, of Tin, where 
Hb = {\v) e TiclEfjo'^-- = 0}' and = {x\v). This 
subspace has dimension 2" — 1 = D — 1, and it is possible 
to choose an orthonormal basis, = 2, . . . D}, for it 

in which each basis element corresponds to a particular 
balanced Boolean function flBI . 

Let us first see how the filtering procedure performs 
when applied to the problem of distinguishing \wk){— 
It/ji)) from the set of the D — 1 orthonormal basis states, 
\vi){= in Tib- We assume their a priori proba- 

bilities to be equal, i.e. rji = rj = (l — t]i)/[D — 1) for 
i = 2, . . . D, where 771 is the a priori probability for \wk)- 
For Ulf = \\wlf = h we obtain h = (2" - l)/22'=-2. 
Then the average overlap, S'fc, between \wk) and the set 
of balanced basis vectors can be written as 

Sk = , (11) 

in terms of fk [llj. The failure probabilities are given 
by Eq. (jH)), using S — Sk and, to good approximation, 
the POVM result holds when l/2'=~2 < jj^^ < 2^-2. 
For example, in the case in which all of the a priori 
probabilities are equal, i.e. 771 = 1/-D, we find that 

QfqM = QsQM = QSQM = (1 + fk)/D. From Fig. [2 
the difference between the POVM and the von Neumann 
measurement is at its largest. To good approximation, 
Qpovm/Qsqm = 4/2''/^, which, for fc > 1, shows that 
the POVM can perform substantially better than the von 
Neumann measurements. 

Now that we know how this procedure performs on the 
basis vectors in Tib, we shall examine its performance on 
any balanced function, i.e. we apply it to the problem of 
distingushing \wk) from the set of all states in Hb that 
correspond to balanced functions. The number of such 
states is = D\/{D/2)\'^ and we again assume their a 
priori probabilities to be equal, 77 = (1 — rii)/N. It can 
be shown (l^ that the average overlap between \wk) and 
the set {\v)} is given by the same expression, Eq. (fTT|l . as 
in the previous case. Therefore, much of what was said 
in the previous paragraph remains valid for this case, as 



well, with one notable difference. The case rji = 1/D 
now does not correspond to equal a priori probability for 
the states but, rather, to a priori weight of the sets that 
is proportional to their dimensionality. In this case it is 
the POVM that performs best. In the case of equal a 
priori probability for all states, 771 = 1/{N +1), we are 
outside of the POVM range of validity and it is the first 
standard quantum measurement (SQMl) that performs 
best. Both the POVM and the SQMl are good methods 
for distinguishing functions in Wk from balanced func- 
tions. Which one is better depends on the a priori prob- 
abilities of the functions. 

Classically, in the worst case, one would have to eval- 
uate a function 2"[(l/2) -I- (1/2'^)] -I- 1 times to determine 
if it is in Wk or if it is an even function. Using quan- 
tum information processing methods, one has a very good 
chance of determining this with only one function eval- 
uation. This shows that Dcutsch-Jozsa-type algorithms 
need not be limited to constant functions; certain kinds 
of biased functions can be discriminated as well. 

Unambiguous state discrimination is a procedure that 
is of fundamental interest in quantum information theory. 
Its only application so far has been to quantum cryptog- 
raphy. The results presented here suggest that related 
methods can also serve as a tool in the development of 
quantum algorithms. 
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